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FOREWORD 


Recent  requirements  for  increased  strength  and  service  life  of  machines 
and  structures  have  been  m^t  by  the  use  of  higher  strength  materials  and  new 
fabrication  and  joining  methods.  Simultaneously,  failures  due  to  fracture  have 
increased  relative  to  those  resulting  from  excessive  deformation.  Frequently 
service  conditions  are  such  that  low  temperature  brittle  fracture,  fatigue 
fracture,  and  high  temperature  creep  rupture  must  be  considered  in  a  single 
system.  National  concern  with  increased  safety,  reliability,  and  cost  has 
focused  atten*'ion  upon  these  problems. 

Methods  are  now  available  to  predict  both  fatigue  crack  initiation  life 
and  crack  propagation  life.  Paradoxically  the  materials  properties  required 
for  long  fatigue  crack  initiation  life  are  incompatible  with  the  requirements 
of  high  fracture  toughness.  Thus,  the  conflicting  design  approaches  and 
requirements  placed  on  the  material  are  confusing  and  often  impossible  to  satisfy. 

Numerous  publications  dealing  with  a  variety  of  fracture  problems  have 
led  to  many  new  and  useful  developments.  However,  the  synthesis  of  the  concepts 
into  methods  for  design,  testing  and  inspection  has  lagged. 

This  program  of  study  is  intended  to  contribute  to  the  integration,  cor¬ 
relation,  and  organization  of  mechanics  and  materials  concepts  and  research 
information  into  a  form  that  will,  permit  enlightened  decisions  to  be  made 
regarding  fracture  control.  Reports  are  in  preparation  in  three  categories: 

1.  Research  reports  designed  to  explore,  study  and  integrate  isolated 
and/or  conflicting  concepts  and  methods  dealing  with  life  prediction, 

2.  Reports  to  introduce  and  summarize  the  state-of-the-art  concepts 
and  methods  in  particular  areas,  and 

3.  Example  problems  and  solutions  intended  to  illustrate  the  use  of 
these  concepts  in  decision  making. 


H.  T.  Corten 
Principal  Investigator 
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SUMMARY 

Fatigue  crack  propagation  is  treated  as  a  succession  of  fatigue  crack 
initiation  events.  The  cyclic  stress  and  strain  experienced  by  each  element 
a.'ead  of  the  crack  tip  is  computed  by  means  of  an  elastic-plastic  analysis  of 
the  crack  tip  region.  The  attendant  cumulative  fatigue  damage  is  computed 
in  terms  of  low  cycle  fatigue  properties  of  the  metal.  This  permits  a  corre¬ 
lation  to  be  made  between  the  fatigue  crack  propagation  resistance  of  a  metal 
and  its  resistance  to  fatigue  as  measured  on  smooth  specimens. 

Over  a  limited  range,  a  log-log  plot  of  the  fatigue  crack  propagation 
rate,  da'dN,  verbis  the  stress  intensity  range,  AK,  is  usually  found  to  be 
nearly  linear  and  may  written  in  the  form:  da/dN  =  A(AK)'\  It  is  found 

in  this  investigat  on  that  the  coefficient,  A,  and  the  exponent,  p,  depend 
upon  the  types  of  assumo^ions  used  in  the  analysis  to  deal  with  the  stress- 
strr.in  singularity  at  the  notch  tip.  If  the  crack  tip  radius  is  assumed  to 
be  blunted  to  the  value  of  the  crack  opening  displacement,  p  is  found  to  be 
equal  to  two  and  A  is  only  a  function  of  the  low  cycle  fatigue  properties  of 
the  material.  If  the  additimal  assumption  is  made  that  each  metal  has  a 
microstructure  size,  p*,  below  which  continuum  mechanics  solutions  are  not 
applicable,  the  coefficient  A  is  found  to  be  a  function  of  AK.  This  causes 
a  significant  deviation  from  the  log-log  linear  relationship  usually  assumed 
to  exist  between  do/dN  and  AK.  The  analysis  indicates  that  at  small  values 
of  AK  the  crack  propagation  rate  rapidly  decreases  until  it  approaches  zero 
at  a  "threshold"  value  of  AK. 

Fatigue  crack  propagation  data  for  eight  steels  from  Barsom's  wc  k  are 
compared  with  predictio.is  from  the  analysis  using  low  cycle  fatigue  properties 
and  values  of  p*  based  on  Gurland's  study  of  the  dimension  of  the  micro- 
structural  elements  responsible  for  strengthening  in  steers.  In  all  cases  the 
estimate  of  c.a/dN  was  within  a  factor  of  two  of  the  measured  values. 
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Strain  hardening  exponent.  Cyclic  value 

Fatigue  crack  propagation  exponent 

Coordinates  along  and  normal  to  crack  length 

Fatigue  crack  propagation  coefficient 

Crack  opening  displacement  at  maximum  load 

Young's  modulus  and  shear  modulus 

Plane  strain  critical  stress  intensity  factor 

CycJes,  Cycles  to  failure.  Transition  fatigue  life 

Reversed  plastic  zone  size 

Mean  free  path  as  defined  by  Gurland 

Effective  crack  tip  radius 

Microstructure  size  related  to  X 

Normal  stress  and  normal  strain 

Fatigue  strength  and  ductility  coefficients 

Monotonic  yield  strength 

Cyclic  uniaxial  yield  strength  and  yield  strain 

Shear  stress  and  shear  strain 

Shear  yield  strength  and  yield  strain 

Stress  intensity  factor  range  (defined  as  AS/iTa  ) 

Threshold  stress  intensity  factor  range 

Nominal  stress  range 

Plastic  strain  range  and  elastic  strain  range 
Normal  stress  range  and  normal  strain  range 
Shear  stress  range  and  shear  strain  range 
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INTRODUCTION 


The  problem  of  plane  strain  fatigue  crack  propagation  has  been  the 
subject  o*^  extensive  study  over  the  last  decade.  In  1961  Paris,  Gomez 
and  Anderson  (1)  suggested  a  correlation  between  the  fatigue  crack  propaga¬ 
tion  race  and  the  range  of  the  applied  stress  intensity  factor.  Experimental 
results  are  now  generally  presented  in  a  form  sirailar  to  Fi)'.  1  which  is  a 
schematic  log-log  plot  of  the  crack  propagation  rate,  da/dN,  as  a  function 
of  the  stress  intensity  factor  range,  AK,  (defined  as  AS  /tt a  )  .  Most 
metals  appear  to  exhibit  a  "threshold"  value  of  AK,  below  which  fatigue 
cracks  do  not  measurably  propagate.  Above  the  threshold  value  there  Ac. 
usually  an  intermediate  range  of  AK  over  which  da/dN  is  nearly  linearly 
related  to  AK  on  a  log- log  basis.  At  higher  ranges,  the  crack  propagation 
rate  increases  rapidly  as  the  maxiirum  stress  intensity  approaches  the  critical 
value  of  K,  .  The  crack  propagation  rate  in  the  intermediate  zone  is  related 
to  the  range  of  stress  intensity  by  the  following  equal  ion: 

si  - 

where  A  and  p  are  experimentally  determined  constants. 

A  variety  of  tests  (2-6)  have  shown  that  for  most  metals  the  exponent 
p  in  Eq.  (l)usually  varies  between  two  and  four,  though  values  as  large  as 
eight  have  occasionally  been  reported.  This  makes  Eq.  (1)  dimensionally 
cumbersome.  Depending  on  the  value  of  the  exponent,  p,  the  constant.  A, 
will  have  different  dimensions. 

Using  continuum  mechanics  variables  in  a  dimensional  analysis,  Liu  (7) 
has  shown  that  the  exponent  p  in  Eq. (1)  should  be  two.  The  experimentally 
observed  deviation  of  p  from  two  can  be  accommodated  within  the  framework 


of  continuum  mechanics  by  permitting  the  coefficient,  A,  in  Eq.  (1)  to  be 


a  function  of  AK.  This  requires  the  introduction  of  a  "length  parameter" 
which  may  be  regarded  as  a  measure  of  the  microstructural  dimension  in  the 
metal  below  which  usual  continuum  mechanics  does  not  apply.  In  tliis  paper 
a  dimensionally  compatible  relationship  is  proposed  between  AK  and  da/dN 
which  has  been  derived  by  analysiing  a  simple  fatigue  crack  propagation 
model  similar  to  that  proposed  by  Liu  and  lino  (8) . 

There  are  practical  as  well  as  theoretical  reasons  for  attempting  to 
relate  fatigue  crack  propagation  to  other  material  properties  via  the  analysis 
of  a  model  such  as  used  here.  A  significant  problem  facing  materials  engi¬ 
neers  is  the  selection  of  materials  to  resist  fatigue  crack  propagation  at 
intermediate  and  low  stresses.  Few  guidelines  are  currently  available  for 
choosing  the  most  crack  resistant  metal  from  a  group  of  candidates  short  of 
conducting  fatigue  crack  propagation  tests.  From  this  point  of  view,  it 
would  be  desirable  to  relate  A  and  p  in  Eq.(l)to  other  known  mechanical 
properties  of  metals. 

Review  of  Fatigue  Crack  Propagacion  Models 

Several  in^^estigators  have  attempted  to  correlate  fatigue  crack  propa¬ 
gation  with  other  mechanical  properties.  For  example,  Tomicins  (9)  has  ana¬ 
lyzed  fatigue  crack  propagation  using  the  plastic  cohesive  stress  approach 
of  Dugdale  (10)  and  concluded  that  fatigue  crack  growth  rate  is  basically 
dependent  upon  the  cyclic  stress-strain  characteristics  of  materials. 

Hickerson  and  Hertzberg  (11)  have  considered  the  correlation  of  fatigue  crack 
growth  resistance  of  materials  with  other  mechanical  properties  and  found  the 
cyclic  strain  hardening  exponent,  cyclic  yield  strength,  and  elastic  modulus 
to  be  important  parameters. 
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A  number  of  continuum  mechanics  models  of  fatigue  crack  propagation 

have  been  analyzed  in  the  past.  For  example,  McCllntock  (12)  proposed  a 

fatigue  crack  propagation  model  on  the  basis  that  fracture  at  the  crack  tip 

occurs  whenever  the  accumulated  average  plastic  strain  reaches  a  critical 

value  in  a  small  sector  of  material  just  ahead  of  the  crack  tip.  The  radius 

of  the  sector  which  McClintock  called  the  structure  size,  was  interpreted 

as  the  smallest  region  to  which  a  macroscopic  fracture  criterion  can  be 

-4 

applied  (he  used  a  value  of  2  x  10  in.).  From  his  damage  analysis, 

McClintock  concluded  that  the  crack  propagation  rate  is  proportional  to  the 
square  of  the  plastic  zone  size. 

Rice  (13)  has  considered  fatigue  crack  growth  based  on  the  plasticity 
model  of  slip  ahead  of  the  crack.  Tracing  the  deformation  history  of  a 
particular  point  after  it  is  encompassed  by  the  plastic  zone  until  the  crack 
tip  advances  to  the  point,  separation  was  assumed  to  occur  when  the  total 
absorbed  hysteresis  energy  equalled  a  postulated  critical  value  per  unit 
area.  He  obtained  a  relation  similar  to  Eq.  (1)  with  p  =  4  but  his  model 
does  not  provide  a  direct  interpretation  of  the  physi(.al  meaning  of  the 
critical  hysteresis  energy. 

Using  Neuber's  rule  (14),  Weiss  (15)  derived  a  crack  propagation  law 
under  the  assumption  that  the  incremental  growth  per  cycle  is  equal  to  the 
distance  over  which  the  strain  ahead  of  the  crack  tip  exceeds  seme  critical 
value.  The  analysis  of  Weiss  included  a  microstructure  length  constant  ,  p*, 
as  introduced  by  Neuber  to  analyze  the  fatigue  behavior  of  sharp  notches. 

Liu  and  lino  (8)  proposed  a  model  in  which  it  is  assumed  that  the 
points  ahead  of  the  crack  tip  constitute  a  set  of  uniaxial  fatigue  specimens. 
The  propagation  of  a  fatigue  crack  is  assumed  to  be  due  to  the  successive 
fatigue  failure  of  each  of  these  imaginary  specimens.  By  postulating  equal 
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crack  advance  per  cycle  and  experimentally  measuring  the  plastic  zone  size 
and  the  strain  distribution  ahead  of  the  crack  tip,  they  derived  an  equation 
similar  to  Eq.  (1)  with  p  =  2. 

Fleck  and  Anderson  (16)  have  analyzed  a  model  similar  to  that  used  by 
Liu  and  lino  (8).  They  used  Rice's  (13)  discrete  surface  of  slip  or  tensile 
yielding  analysis  (for  a  perfectly  plastic  material)  to  estimate  the  plastic 
zone  size  ahead  of  the  crack  tip.  It  was  postulated  that  a  fatigue  element 
fails  when  the  crack  has  progressed  half  way  through  the  element.  By  suitably 
choosing  the  uniaxial  material  properties  and  the  element  size,  Fleck  and 
Anderson  (17)  obtained  correlation  between  experimental  data  and  their  anal¬ 
ysis  for  various  metals. 

Scope 

The  present  paper  adopts  Liu's  model,  but  instead  of  experimentally 
determining  the  plastic  zone  size  and  the  strain  distribution  ahead  of  the 
crack,  the  analytical  results  obtained  by  Rice  (18)  are  employed.  In  its 
simplest  form,  the  analysis  results  in  the  exponent,  p,  of  Eq.  (1)  being  equal 
to  two.  The  coefficient.  A,  can  be  calculated  from  the  usual  cyclic  stress- 
strain  and  fatigue  properties  of  uniaxial  specimens  of  a  given  metal. 

Further,  by  introducing  a  microstructure  size,  p* ,  into  the  analysis  it  is 
shown  that  a  nonlinear  relationship  between  log  AK  and  log  da/dN  may  be 
obtained  which  quantitatively  agrees  with  the  trends  in  experimental  data. 

The  nonlinearity  is  achieved  by  keeping  the  power,  p,  equal  to  two  while 
the  coefficient.  A,  becomes  a  function  of  AK.  Thus,  the  dimensionality  of 
Eq.  (1)  is  maintained. 


DESCRIPTION  OF  THE  FATIGUE  CRACK  PROPAGATION  MODEL 


Figure  2  shows  the  model  to  be  analyzed  and  a  schematic  representation 
of  the  stress  and  strain  history  experienced  by  an  element  ahead  of  a  pro¬ 
pagating  crack  in  a  plate  subjected  to  repeated  tension  so  that  AK  is  constant. 
Figure  3a  depicts  an  essentially  undamaged  element  ab  at  the  reversed 
plastic  zone  boundary.  Figure  3b  shows  the  same  element  with  the  crack  at 
a  distance  x  from  it.  The  damage  experienced  by  element  ab  in  one  cycle 
at  this  location  will  be  the  same  as  experienced  by  element  cd  in  Fig.  3a 
(assuming  AK  is  constant).  Thus,  the  cumulative  damage  required  to  cause 
fatigue  failure  of  an  element  is  equivalent  to  the  sum  of  the  damage  experi¬ 
enced  during  one  cycle  by  each  element  lying  within  the  reversed  plastic  zone. 
Let  N^(x)  be  the  fatigue  life  of  a  specimen  of  the  material  subjected  to 
constant  cyclic  strain  equal  to  that  which  is  present  at  a  distance  x  from 
the  crack  tip.  Using  linear  fatigue  damage  summation,  failure  of  the  element 
will  occur  when 

R 

P 

^  ^  N^  (2) 

X  =  0 


Assuming  that  AK  is  constant  the  crack  advances  in  each  cycle  a  con¬ 
stant  amount.  Ax  =  da/dN.  Mulitplying  the  left  side  of  Eq.(2)  by  da/dN  and 
the  right  side  by  Ax  and  approximating  the  sum  by  an  integral: 


da 

dN 


N^(x) 


dx 


(3) 


Approximating  the  sum  by  an  integral  is  equivalent  to  approximating  the  area 
under  the  step-like  curve  by  the  area  under  the  smooth  dotted  curve  in  Fig.  3c. 


MECHANICS  AND  TATIGUE 


ANALYSIS  OF  THE  CRACK  TIP  REGION 

Rice  (13,  18)  considered  a  stationary  crack  loaded  in  antiplane  shear 
(mode  III)  under  small  scale  yielding  in  a  material  with  the  following  shear 
stress-strain  law: 


t“Gy  forY<Y,T<T 

—  0—0 

T  =  T  (y/y  )^  for  Y  ^  Y  »  ^ 

o  o  —  o  —  o 


(4) 


Assuming  that  the  plastic  strain  components  at  each  point  within  the 
plastic  zone  remain  proportional  to  each  other,  Rice  derived  the  following 
expressions  for  cyclic  stress  and  strain  ranges: 


At  ■  2t 

o 


AK 


III 


4(1  +  n)  •JT  T  X 
o 


1  +  n 


Ay  *  2y, 


AK 


III 


4(1  +  n)  IT  T  X 
o 


1  +  n 


(5) 


A  similar  analysis  for  tensile  loading  (mode  I)  is  not  available. 

McClintock  (12)  has  discussed  the  analogy  between  mode  III  and  mode  I 
for  the  case  where  displacements  parallel  to  the  crack  are  small  compared  to 
those  normal  to  the  crack  surface.  In  the  present  study  the  above  analogy 
will  be  used  to  approximate  the  stress  and  strain  ranges  ahead  of  a  tensile 
crack  in  a  material  obeying  the  following  cyclic  stress-strain  law: 
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(6) 


In  terms  of  the  product  of  stress  and  strain  range,  Eq,  (5)  reduces  to 


Ao  Ae 


AK‘ 


(1  +  n' )  IT  Ex 


(7) 


The  reversed  plastic  zone  size,  R^,  and  the  crack  opening  displacement 
at  maximum  load,  COD,  are  given  by  the  following: 


R 

P 


1  _ 

4(1  +  n') 


COD 


2  e'  AK 

_ 1 _ 


TT  O' 


(8) 


(9) 


Even  though  the  above  equations  are  derived  for  a  stationary  crack, 
they  will  be  assumed  to  hold  for  a  propagating  crack.  The  use  of  the  fully 
stabilized  cyclic  stress-strain  law  may  be  justified,  provided  the  crack 
propagation  rate  is  slow  enough  that  each  point  within  the  plastic  zone 
spends  sufficient  number  of  cycles  at  each  strain  range  to  attain  stable 
stress-strain  behavior. 

In  order  to  compute  the  fatigue  life  of  an  element  ahead  of  the  crack 
tip  subjected  to  a  known  stress-strain  range,  the  following  power  functions 
will  be  used: 


AO  -  2o^  (2N^)‘’ 

(10a) 

Ae  =  2eJ  (2N^)‘^ 

(10b) 

Note  that  in  Eq.  (10b) the  fatigue  life  has  been  assumed  to  be  related 
to  the  total  strain  range  by  a  power  function.  In  most  materials,  however. 


the  total  strain  range  is  related  to  the  fatigue  life  by  the  following 
equation  (Fig.  4); 


f  V 

Ae  =  Ae  +  Ae  »  2e'  (2N^)‘^  +  — ^  (2Nj  (11a) 

p  e  t  t  E  i: 

Equation  (11a),  in  general,  cannot  be  solved  for  the  fatigue  life 
explicitly  and  consequently  a  numerical  method  has  to  be  adopted.  On  the 
other  hand  a  closed  form  solution  to  the  present  problem  may  be  obtained  if 
the  following  approximations  are  made.  The  majority  of  the  damage  experienced 
by  a  fatigue  element  ahead  of  the  crack  occurs  well  inside  the  reversed 
plastic  zone  (17)  where  the  plastic  strain  range  is  much  larger  than  the 
elastic  strain  range.  Consequently,  in  the  low  cycle  (high  strain)  regime 
of  fatigue,  the  fatigue  life  for  a  given  strain  range  will  lie  between  the 
following  two  power  functions  with  slightly  different  exponents,  c  and  c’, 
as  defined  in  Fig.  4: 


Ae  Ae  =  2e!. 
P  f 

(2Np‘' 

(11b) 

Ae  2e^ 

(2Np‘'' 

(11c) 

The  exponent  c  is  the  usual  fatigue  ductility  exponent  or  the  slope 
of  the  log  ACp  versus  log  2N^  plot.  A  straight  line  on  the  log-log  plot 
connecting  the  intercept  at  one  reversal  (the  fatigue  ductility  coefficient, 
ep  and  the  strain  range  at  the  transition  fatigue  life  will  have  a  slope  of 
c ' . 

Note  that  Eq. (11b)  always  underestimate r  and  Eq. (11c)  always  overestimates 
the  fatigue  life  in  the  relevant  range  of  the  curve  in  Fig.  4.  Thus  the  crack 
propagation  analysis  based  on  Eqs.  (11b)  and  (11c)  should  bracket  the  results 
that  would  be  obtained  if  the  unabridged  Eq.  (11a)  were  used  to  compute 
fatigue  lives. 


Multiplying  Eqs.  (10a)  and  (10b)  and  equating  the  result  to  Eq.  (7) 
gives  the  damage  sustained  in  one  cycle  by  an  element  at  a  distance  x 
from  the  crack  tip. 


-  =  2 
"f 


AK 


4(1  +  n')  Tt  E  a[T[ 


-1  1 

b+c  b+c 

(x) 


(12y 


A  disturbing  aspect  of  Eq.  (12)  is  that  the  integral  in  Eq.  (3)  is 
divergent  for  usual  values  of  b  and  c.  This  singularity  can  be  avoided 
by  introducing  crack  tip  blunting.  iX  any  instant  the  crack  tip  will  have 
a  finite  radius  that  is  approximately  equal  to  half  of  the  crack  opening 
displacement.  To  treat  this  problem  rigorously,  a  full  nonlinear  analysis 
is  needed  which  unfortunately  is  not  a.-ailable.  In  the  absence  of  such  an 
analysis  Rice's  results  are  modified  by  translating  the  origin  of  the  axes 
into  the  crack  a  distance  of  COD/2.  In  other  words,  x  in  Eq.  (12)  is 
replaced  by  x  +  COD/2.  Substituf!  ig  the  modified  form  of  Eq.  (12)  into 
Eq.  (3)  and  using  Eqs,  (8)  and  (9)  gives: 


da 

dN 


-2(b+c) 

b+c+1 


4(1  +  n')  Oj 


-1 

b+c 


1  -  [4(1  +  n')ej^] 


b+c+1 1 

e'  AK^ 

1  y 

b+c  j 

J 

,2 

Ti  a 

y 

(13) 


Since  all  terms  in  Eq.  (13)  except  da/dN  and  AK  are  constants  for 
a  given  metal,  it  may  be  written  in  a  form  similar  to  Eq.  1  with  p  =  2  as 
follows: 

^  =  A(AK)^  (13a) 

dN 

The  coefficient  A,  is  a  function  of  the  cyclic  stress-strain  and  fatigue 
properties  of  the  metal.  In  the  derivation  of  Eq.  (13),  o^,  and  n'  have 
been  used  to  characterize  the  cyclic  deformation  resistance  and  o^,  e^,  b 
and  c,  the  fatigue  resistance  of  the  metal. 
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INFLUENCE  OF  MATERIAL  PROPERTIES  ON  THE  COEFFICIENT  OF  EQ.(13) 

Equation  (13)  requires  the  exponent  of  AK  in  Eq.  (1)  to  be  two  and 
the  coefficient  A  is  expressed  as  a  function  of  the  cyclic  stress-strain 
and  fatigue  properties  of  materials.  For  any  fixed  AK,  the  A's  are  pro¬ 
portional  to  the  rate  of  fatigue  crack  propagation  for  different  metals. 

Thus  the  quantitative  influence  of  these  properties  on  fatigue  crack  propa¬ 
gation  resistance  may  be  evaluated  by  calculating  their  effect  on  tl  e 
coefficient  A. 

For  making  such  an  evaluation,  reference  values  for  t',  o'  o',  n' , 

f  f  y 

b,  c  and  E  shown  in  Fig.  5  were  chosen  and  each  property  in  turn  was  varied 
keeping  the  others  fixed.* 

The  effects  of  the  various  properties  on  the  coefficient  A  are  shown 
in  dimensionless  form  in  Fig.  5.  Though  the  variation  in  one  material  property 
without  accomnanying  variation  in  the  others  is  highly  unrealistic,  these 
plots  still  serve  to  indicate  how  sensitive  the  coefficient  A  is  to  small 
changes  in  material  properties.  This  helps  separate  those  properties  that 
need  to  be  known  accurately  from  those  that  need  not  be. 

It  is  seen  in  Fig.  5  that  the  coefficient  A  is  more  sensitive  to 
relative  changes  in  fatigue  strength  coefficient  (o^,  fatigue  ductility 
coefficient  and  exponent  (e^  and  c)  and  elastic  modulus  (E)  than  to 
changes  in  other  material  properties.  These  four  properties,  therefore,  have 
to  be  known  more  accurately  than  the  rest  for  proper  evaluation  of  fatigue 
crack  resistance. 

It  is  interesting  to  note  that  according  to  the  present  analysis,  the 
cyclic  yield  strength  (o^)  and  the  cyclic  strain  hardening  exponent  (n') 
appear  to  be  of  relatively  minor  importance  in  determining  fatigue  crack 
propagation  resistance. 

* 

The  reference  values  in  Fig.  5  were  chosen  to  be  representative  nt  a 
slightly  hardened  structural  steel. 


EFFECT  OF  MICRO STRUCTURE  SIZE 


A  microstructure  size,  p*,  may  be  incorporated  in  the  analysis  as  follows. 
Assume  there  is  a  distance  p*  ahead  of  the  crack  tip  that  is  a  function  of 
the  microstructure  of  the  material  within  which  the  usual  continuum  mechanics 
apprr/ach  is  not  meaningful.  In  other  words,  continuum  analysis  and  bulk 
fatigue  properties  are  applicable  only  from  a  distance  p*  ahead  of  the  crack 
tip,  and  Eq.(3)  should  be  modified  by  changing  the  lower  limit  of  the  integral 
from  0  to  p*.  The  results  will  then  be  applicable  provided  the  reversed 
plastic  zone  size  R^  is  large  compared  to  p*.  Incorporating  these  modifi¬ 
cations  into  the  analysis  Eq.  (13)  is  altered  as  follows: 


da  _  -2 (b+c) 

r  'j'  1 

y 

-1 

b+c 

4 

1+  M 

COD 

b+c+1 

b+c 

4(l+n* 

dN  b+c+1 

4(l+n')aje^ 

\ 

1 

I 


b+c+1 

b+c 

> 


e'AK^ 

JL _ 

y 


(14) 


The  condition  R^  >>  p*  implies 

AK  »  2(1  +  n')^^^  o^/if^  (15) 

Notice  that  for  large  enough  load  levels,  COD  will  be  large  compared  to  p* 
and  Eq.  (14)  becomes  the  same  as  Eq.  (13).  For  small  loads  where  COD  is 
of  the  order  of  p*,  Eq.  (14)  will  give  significant  reduction  in  da/dN  from 
Eq.  (13).  Thus,  in  general,  the  log  da/dN  versus  log  AK  plot  will  be 
nonlinear  and  will  asymptotically  approach  the  linear  relationship  of  Eq.  (13) 
for  large  values  of  AK.  At  small  AK's  the  plastic  zone  size  will  approach 
the  value  of  p*  which,  according  to  this  analysis,  will  arrest  crack 
growth.  Thus,  the  limiting  value  of  AK  from  Eq.  (15)  may  be  regarded  as 
an  estimate  of  the  threshold  ’'alue,  4K^j^. 

Figure  6  contains  plots  of  AK  versus  da/dN  for  a  typical  steel 
(same  properties  as  used  as  reference  in  Fig.  5)  and  for  various  values  of 


12 


p*.  It  may  be  observec  that  the  effect  of  p*  is  significant  at  low  values 
of  AK  and  tends  to  bf  less  significant  at  higher  values  of  AK.  The  fatigue 
crack  growth  resistance  of  any  metal  may  then  be  thought  of  as  consisting  if 
two  parts.  Firstly,  the  fatigue  crack  growth  resistance  is  controlled  by 
the  bulk  low  cycle  fatigue  properties  and  is  given  by  Eq.  (13).  Secondly, 
an  increase  in  resistance  beyond  that  predicted  by  Eq.(13)  is  achieved  by 
introducing  the  microstructure  size,  p*.  It  is  the  latter  effect  that  makes 
it  necessary  to  select  an  exponent  p  of  Eq.  (1)  that  is  greater  than  two 
in  order  to  fit  the  power  function  to  experimental  data  over  a  limited  range 
of  AK. 

The  physical  interpretation  of  tlie  microstructure  size,  p*,  is  not 
clear.  Some  possibilities  would  include  the  distance  between  slip  bands, 
grain  size,  carbide  spacing,  mean  free  ferrite  path,  and  so  on. 

The  effects  of  size  or  scale  of  the  microstructure  on  the  yield  strength 
of  steels  have  been  studied  in  some  detail.  Recently,  Gurland  (19)  has  shown 
that  in  many  steels  of  various  microstructures,  the  yield  strength  can  be 
expressed  as  a  function  of  a  microstructu*'"!  parameter,  X,  in  a  form  similar 
to  the  Hall-Petch  equation,  i.e.: 


0 

y 


0  +  k 

o  y 


(16) 


0  *  Friction  stress 


k  “  Constant 

y 


fcrain  size  in  low  carbon  mild  steels 


X  *<Mean  free  path  between  cementite  par- 

[ticles  for  spheroidic  and  tempered  steels 


Mean  cell  diameter  in  deformed  pearlitic 
(^steels 


When  X  is  interpreted  as  above,  data  for  a  variety  of  steels  correlate 
well  with  a  single  Hall-Petch  equation.  For  the  purposes  of  the  present 
paper  p*  for  steels  will  be  approximated  by  X. 


/ 

/ 

13  / 

/ 

f 

COMPARISON  OF  EQ.  (14)  WIt/ BARSOM'S  DATA  ON  STEEL 

/ 

/ 

Published  crack  propagation  /data  are  seldom  accompanied  by  cyclic 
or  low  cycle  fatigue  properties  of  the  metals  studied.  To  compare  crack 
propagation  r'.ata  from  the  literature  with  Eq.  (14)  ,  therefore  requires 
that  the  arpropriate  mechanical  properties  first  be  estimated.  Birsom's 
data  (6)  on  fatigue  crack  propagation  in  steels  were  chosen  for  comparison 
with  Eq.  (14).  The  fatigue  strength  and  ductility  coefficients  were  approx¬ 
imated  from  the  reported  values  of  ultimate  tensile  strengtli  and  true  frac¬ 
ture  ductility  (20).  The  fatigue  ductility  exponent,  c,  was  approximated 
as  -0.60  for  the  four  soft  steels  and  as  -0.65  for  the  harder  steels. 

The  various  properties  assumed  for  computation  are  shown  in  Table  X. 

As  ms^tioned  before,  p*  has  been  approximated  by  X  as  given  in 
Eq.  (16).  Curland's  (19)  reported  values  for  =  10.5  ksi  and  k^  =  0.564 
ksi  /in  were  used  in  Eq.  (16) . 

Figures  7  through  14  show  the  predicted  values  of  crack  propagation 
rate  according  to  Eq.  fl4)  together  with  Barsom's  test  data.  The  predicted 
values  have  been  based  on  two  values  of  c  in  Eq.  (10b)  as  discussed  earlier. 
For  all  practical  purposes,  Barsom's  test  data  fall  within  the  band  of  pre¬ 
dicted  values  in  all  cases.  Notice  that  the  value  of  p*  for  the  high 
yield  steels  are  on  the  average  a  factor  of  ten  smaller  than  those  for  the 
low  yield  steels.  Over  the  range  of  AK  investigated  by  Barsom  for  the 
high  yield  strength  steels,  the  effect  of  p*  is  small.  On  the  other  hand, 
for  the  low  yield  strength  steels  the  introduction  of  p*  gives  .-xcellent 
correlation  with  both  the  magnitude  of  the  crack  propagation  rate  and  the 
slope  of  the  log  AK  -  log  da/dN  plots. 
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CONCLUDING  REMARKS 

Ignoring  the  effect  of  micro structure  size,  Eq.  (13)  gives  the  maximum 
potential  rate  of  fatigue  crack  propagatioi.  for  any  given  range  of  stress 
intensity  factor,  and  the  slope  of  the  log  AK  -  log  da/dN  plot  is  two.  This 
maximum  potential  fatigue  crack  propagation  rate  is  controlled  by  the  bulk 
mechanical  properties  o^  the  metal,  particularly  the  fatigue  ductility 
coefficient  and  exponent,  the  fatigue  strength  coefficient,  and  the  '.lastic 
modulus.  For  this  kind  of  fatigue  crack  resistance  the  yield  strength  is 
theoretically  of  little  importance  (Fig.  5). 

By  introducing  a  microstructure  size  into  the  analysis  the  crack  propa¬ 
gation  rate  for  any  given  stress  intensity  factor  range  is  reduced,  especially 
at  intermediate  to  low  values  of  stress  intensity  factor.  The  fncrease  in 
fatigue  crack  propagation  resistance  is  larger  the  larger  the  microstructure 
size  (Fig.  6).  In  steels,  a  Hall-Petch  type  equation  may  be  used  for  com¬ 
puting  the  relevant  microstructure  size  from  the  yield  strength.  Since  the 
microstructure  size  varies  inversely  as  the  yield  strength,  this  kind  o 
fatigue  crack  propagation  resistance  is  indirectly  dependent  on  the  yield 
•"trength  of  steels. 
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Fig.  2  SCHEMATIC  OF  STRESS- STRAIN  HISTORY 
OF  A  FATIGUE  ELEMENT  AHEAD  OF  THE 
CRACK  TIP 


Reversed  Plastic  Zone 


(a)  ELEMENT  AT  EDGE  OF  REVERSED  PLASTIC  ZONE 

Reversed  Plastic  Zone 


(b)  ELEMENT  AT  DISTANCE  x  FROM  CRACK  TIP 


Distance  from  Crack  Tip 
(c)  DAMAGE  HISTORY  OF  AN  ELEMENT 


Fig.  3  TYPICAL  DAMAGE  HISTORY  OF  AN  ELEMENT 


log(2NJ 


Fig.  4  APPROXIMATING  TOTAL  STRAIN  AMPLITUDE 
VERSUS  LIFE  PLOTS  BY  POWER  FUNCTIONS 
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Fig.  7  COMPARISON  OF  EXPERIMENTAL  DATA(6)  WITH 
THEORETICAL  PREDICTION  FOR  HY-80  STEEL 
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Fig.  8  COMPARISON  OF  EXPERIMENTAL  DATA(6)  WITH 
THEORETICAL  PREDICTION  FOR  HY-130  STEEL 
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Fig.  9  COMPARISON  OF  EXPERIMENTAL  DATA(6)  WITH 
THEORETICAL  PREDICTION  FOR  lONi-Cr-Mo-Co  STEEL 
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Fig,  10  COMPARISON  OF  EXPERIMENTAL  DATA(6)  WITH 
THEORETICAL  PREDICTION  FOR  l2Ni-5Cr-3Co  STEEL 
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Fig.  13  COMPARISON  OF  EXPERIMENTAL  DATA(6)  WITH 
THEORETICAL  PREDICTION  FOR  A302-B  STEEL 


Fig.  14  COMPARISON  OF  EXPERIMENTAL  DATA(6)  WITH 
THEORETICAL  PREDICTION  FOR  A537-A  STEEL 


